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The STAR collaboration reports the first observation of exclusive p° photo- production, AuAu^ 
AuAup , and p° production accompanied by mutual nuclear Coulomb excitation, AuAu —* 
Au*Au*p°, in ultra-peripheral heavy-ion collisions. The p° have low transverse momenta, con- 
sistent with coherent coupling to both nuclei. The cross sections at ^/snn = 130 GeV agree with 
theoretical predictions treating p° production and Coulomb excitation as independent processes. 
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In ultra-peripheral heavy-ion collisions the two nuclei 
geometrically 'miss' each other and no hadronic nucleon- 
nucleon collisions occur. At impact parameters b sig- 
nificantly larger than twice the nuclear radius Ra, the 
nuclei interact by photon exchange and photon-photon 
or photon-Pomeron collisions 13]. Examples are nuclear 
Coulomb excitation, electron-positron pair and meson 
production, and vector meson production. The exchange 
bosons can couple coherently to the nuclei, yielding large 
cross sections. Coherence restricts the final states to low 
transverse momenta, a distinctive experimental signa- 
ture. The STAR collaboration reports the first observa- 
tion of coherent exclusive p° photo-production, AuAu—> 
AuAup , and coherent p° production accompanied by 
mutual nuclear excitation, AuAu — ► Au*Au*p°. Ultra- 
peripheral heavy-ion collisions are a new laboratory for 
diffractive interactions, complementary to fixed-target p° 
photo-production on complex nuclei Q. 

Exclusive p° meson production, AuAu^AuAup (c.f. 
Fig. [l^,), can be described by the Weizsacker- Williams 
approach || to the photon flux and the vector meson 
dominance model [Q. A photon emitted by one nucleus 
fluctuates to a virtual p° meson, which scatters elasti- 
cally from the other nucleus. The gold nuclei are not 
disrupted, and the final state consists solely of the two 
nuclei and the vector meson decay products |5). In the 



rest frame of the target nucleus, mid-rapidity p° produc- 
tion at RHIC corresponds to a photon energy of 50 GeV 
and a photon-nucleon center-of-mass energy of 10 GeV. 
At this energy, Pomeron (V) exchange dominates over 
meson exchange, as indicated by the rise of the p° pro- 
duction cross section with increasing energy in lepton- 
nucleon scattering M . In addition to coherent p° produc- 
tion, the exchange of virtual photons may excite the nu- 
clei. These processes are assumed to factorize for heavy- 
ion collisions, which is justified by the similar case of 
two-photon interactions in relativistic ion collisions ac- 
companied by nuclear breakup, where it was shown that 
the non-factorisable diagrams are small §. The process 
AuAu-^ Au*Au*p° is shown in Fig. [l]b. In lowest order, 
mutual nuclear excitation of heavy ions occurs by the ex- 
change of two photons ||, Because of the Coulomb 
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FIG. 1: Diagram for (a) exclusive p° production in ultra- 
peripheral heavy ion collisions, and (b) p° production with 
nuclear excitation. The dashed lines indicate factorization. 
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barrier for the emission of charged particles, nearly all 
nuclear decays following photon absorption include neu- 
tron emission fic(| . 

The photon and Pomeron can couple coherently to 
the gold nuclei. The wavelength X lt -p > 2Ra leads to 
coherence conditions: a low transverse momentum of 
p T < Trh/R A (~ 90 MeV/c for gold with R A ~ 7 fm), 
and a maximum longitudinal momentum of p\\ < nhj/ Ra 
(~ 6 GeV/c at 7 = 70), where 7 is the Lorentz boost of 
the nucleus. The p° production cross sections are large. 
The photon flux is proportional to the square of the nu- 
clear charge Z 2 || , and the forward cross section for elas- 
tic p°A scattering da pA / dt\t=o scales as A 4 / 3 for surface 
coupling and A 2 in the bulk limit. At a center-of-mass 
energy of -^/SaT/v = 130 GeV per nucleon-nucleon pair, a 
total p° cross section, regardless of nuclear excitation, 
a(AuAu — > Au^Au^p ) = 350 mb is predicted from a 
Glauber extrapolation of ^p~^ p°p data ||. Calculations 
for coherent p° production with nuclear excitation as- 
sume that both processes are independent, sharing only 
a common impact parameter ^ || . 

In the year 2000, the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory collided gold 
nuclei at ^s^n = 130 GeV. In the Solenoidal Tracker at 
RHIC (STAR) pH , charged particles are reconstructed 
with a cylindrical time projection chamber (TPC) jl^ ] 
operated in a 0.25 T solenoidal magnetic field. A cen- 
tral trigger barrel (CTB) of 240 scintillator slats sur- 
rounds the TPC. Two zero degree hadron calorimeters 
(ZDCs) at ± 18 m from the interaction point are sensi- 
tive to the neutral remnants of nuclear break-up, with 
98 ±2% acceptance for neutrons from nuclear break-up 
through Coulomb excitation || |l3| ]. 

Exclusive p° production has a distinctive signature: 
the 7r + 7r~ from the p° decay in an otherwise 'empty' de- 
tector. The tracks are approximately back-to-back in the 
transverse plane due to the small pr of the pair. The gold 
nuclei remain undetected within the beam. 

Two data sets are used in this analysis. For AuAu—i- 
AuAup , about 30,000 events were collected using a low- 
multiplicity 'topology' trigger. The CTB was divided in 
four azimuthal quadrants. Single hits were required in 
the opposite side quadrants; the top and bottom quad- 
rants acted as vetoes to suppress cosmic rays. A fast 
on-line reconstruction ]l4| removed events without re- 
constructible tracks from the data stream. To study 
AuAu — > Au*Au*p°, a data set of about 800,000 'min- 
imum bias' events, which required coincident detection 
of neutrons in both ZDCs as a trigger, is used. 

Events are selected with exactly two oppositely 
charged tracks forming a common vertex within the in- 
teraction region. The p° candidates are accepted within 
a rapidity range \y p \ <1. A systematic uncertainty of 5% 
is assigned to the number of p° candidates by varying the 
event selection criteria. The specific energy loss dE/dx 
in the TPC shows that the event sample is dominated by 
pion pairs. Without the ZDC requirement in the topol- 
ogy trigger, cosmic rays are a major background. They 
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FIG. 2: The pr spectra of pion pairs for the 2-track events 
selected by (a) the topology trigger (On, On) and (b) the mini- 
mum bias trigger (xn,xn). Points are oppositely charged pairs, 
and the shaded histograms are the normalized like-sign com- 
binatorial background. The open histograms are simulated 
p° superimposed onto the background. 



are removed by requiring that the two pion tracks have 
an opening angle of less than 3 radians. Using the en- 
ergy deposits in the ZDCs, we select events with at least 
one neutron (xn,xn), exactly one neutron (In, In), or no 
neutrons (On, On) in each ZDC, and events with at least 
one neutron in exactly one ZDC (xn,0n); the latter two 
occur only in the topology trigger. A 10% uncertainty 
arises from the selection of single neutron signals. 

The uncorrected transverse momentum spectra of pion 
pairs for the two-track event samples of the topology 
trigger (On, On) and the minimum bias trigger (xn,xn) 
are shown in Fig. |^. Both spectra are peaked at px ~ 
50 MeV/c, as expected for coherent coupling. A back- 
ground model from like-sign combination pairs, normal- 
ized to the signal at pr > 200 MeV/c, is not peaked. For 
comparison, the px spectra from Monte Carlo simula- 
tions H discussed below are shown. They are normal- 
ized to the p° signal at pr < 150 MeV/c and added to 
the background. The M WT invariant mass spectra (c.f. 
Fig. |J) for both event samples are peaked around the p° 
mass. We find 131 ± 14 (On, On) and 656 ± 36 (xn,xn) 
events at px < 150 MeV/c, which we define as coherent 
p° candidates. 

The data contain combinatorial background contri- 
butions from grazing nuclear collisions and incoherent 
photon-nucleon interactions, which are statistically sub- 
tracted. Incoherent p° production, where a photon in- 
teracts with a single nucleon, yields high pt p°s, which 
are suppressed by the low pair px requirement; the 
remaining small contribution is indistinguishable from 
the coherent process. A coherently produced back- 
ground arises from the mis-identified two-photon process 
AuAu Au^ Au^*H + 1~ . It contributes mainly at low 
invariant mass M„ < 0.5 GeV/c 2 . Electrons with mo- 
menta p < 140 MeV/c can be identified by their energy 
loss dE/dx. About 30 e + e~ pairs, peaked at low pair 
p T ^20 MeV/c, were detected in the minimum bias data 
sample |y| . They are extrapolated to the full phase space 
using a Monte-Carlo simulation that describes e + e" pair 
production by lowest order perturbation theory Jl6[ . 
Electron-positron pairs contribute 4 ± 1% to the signal 
at p T < 150 MeV/c and M p ± 0.3 GeV/c. For a given 
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FIG. 3: Rapidity distribution (a) of p° candidates (xn,xn) 
for the minimum bias data (points) compared to the normal- 
ized reconstructed (shaded histogram) and generated (open 
histogram) events from the Monte Carlo simulation. The dif- 
ferential cross section (b) dafaAu—* pAu) /dt for the same data 
set; the line indicates the exponential fit. 



Mu, mucins have lower momenta than the correspond- 
ing electrons and are less likely to be detected. Their 
< 2% contribution to the coherent signal, as well as the 
contribution from u> decays are neglected. 

The acceptance and reconstruction efficiency were 
studied using a Monte Carlo event generator that re- 
produces the expected kinematic and angular distribu- 
tions for p production with and without nuclear excita- 
tion || [It]], coupled with a full detector simulation. The 
p° decay angle distribution is consistent with s-channel 
helicity conservation. The p° production angles are not 
reconstructed since the AuAu scattering plane can not 
be determined. The efficiencies are almost independent 
of pt and the reconstructed invariant mass M^. For the 
minimum bias trigger, 42±5% of all p° within \y p \ < 1 are 
reconstructed. The topology trigger vetoes the top and 
bottom of the TPC, reducing the geometrical acceptance. 
Pions with p T < 100 MeV/c do not reach the CTB, ef- 
fectively excluding pairs with M vv < 500 MeV/c 2 . Only 
7±1% of all p° with \y p \ < 1 are reconstructed in the 
topology trigger. The pr resolution is 9 MeV/c. The 
M-n-Tr and rapidity resolutions are 11 MeV/c 2 and 0.01. 

The rapidity distribution for p° candidates (xn,xn) 
from the minimum bias data is shown in Fig. |3|a). It 
is well described by the reconstructed events from a sim- 
ulation, which includes nuclear excitation [j5|. The gener- 
ated rapidity distribution is also shown. The acceptance 
is small for \y p \ > 1, so this region is excluded from the 
analysis. Cross sections are extrapolated from \y p \ < 1 to 
the full 47r acceptance by a^/a^ y ^ = 1.9 for p° produc- 
tion with nuclear break-up, and a^/a^ y ^ = 2.7 for p° 
production without nuclear break-up. A 15% uncertainty 
in the extrapolations is estimated by varying the Monte 
Carlo parameters. Event rapidity and photon energy are 
related by y = (1/2) In (2E 7 /M p ). But, the average pho- 
ton energy per rapidity bin (E 7 ) ~ 50 GeV is constant, 
when taking the ambiguity of photon emitter and scat- 
tering target into account. 

The minimum bias data sample has an integrated lu- 
minosity of L = 59 mb _1 . The luminosity was mea- 
sured by counting events containing more than 5 nega- 
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FIG. 4: The da(AuAu — > Au* Au*p)/dM^ spectrum for 
2-track (xn,xn) events with pair-pr < 150 MeV/c in the min- 
imum bias data. The shaded histogram is the combinato- 
rial background, and the hatched histogram contains an ad- 
ditional contribution from coherent e + e~ pairs. The fits cor- 
respond to Eq. |j| the sum (solid) of a Breit-Wigner, a mass- 
independent contribution from direct tt + tt~ production and 
their interference (all dashed) , and a second order polynomial 
for the residual background (dash-dotted). 



tively charged hadrons with pt > 100 MeV/c and pseudo- 
rapidity I77I < 0.5. These events represent 79% of the 
hadronic cross section [|l8) . We assume a total gold-gold 
hadronic cross section of 7.2 b Q; its uncertainty domi- 
nates the 10% systematic uncertainty of L. 

The differential cross section da(jAu — > pAu)/dt ~ 
da{^jAu^rpAu)jdp\ for the (xn,xn) events is shown in 
Fig. ||b). Here, the combinatorial background is sub- 
tracted. The photon flux is determined by integration 
of the photon-spectrum of a relativistic nucleus over 
the impact parameter space [^|. For p° production in 
ultra-peripheral collisions, da j dt reflects not only the nu- 
clear form factor, but also the photon px distribution 
and the interference of production amplitudes from both 
gold nuclei. The interference arises since both nuclei 
can be either the photon source or the scattering tar- 
get (jl9). A detailed study of this effect is beyond the 
scope of this paper and the available statistics. From a 
fit to da pAu /dt cc e~ ht we obtain a forward cross section 
dcrP A /dt\ t=0 = 965 ± 140±230 mb/GeV 2 and an approxi- 
mate gold radius of Rau = V4h = 7.5±2 fm, comparable 
to previous results 0. 

The da(AuAu — > Au*Au*p)/dM 1T7r invariant mass 
spectrum for the (xn,xn) events with a pair px < 
150 MeV/c is shown in Fig. [|; the (On, On) events have a 
similar da/dM^ spectrum. Three different parametcri- 
zations are applied: 



da/dM^n = 
da j dM-n-n — 



f p ■ BW(M nw ) + fj ■ I{M nw ) + f p , (1) 

2 



.4 



K7T J\dpY' p 



M 2 — M 2 



iMpTp 



B 



+ f P , (2) 

da/dM^ = f P -BW{M^)-{m p /M^) n + f p . (3) 



Eq. n] is a relativistic Breit-Wigner, BW = M^^M p T p / 
pf;-M 2 J 2 + M 2 r?], for p° production plus a Soding 
interference term M, J(M„) = (A/ 2 — M 2 7r )/[(M 2 — 
A/ 2 ^) 2 +M 2 r 2 1, EqTgis a modified Soding parametriza- 
tion |2l[] , and Eq. || is a phenomenological Ross-Stodolsky 
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Eq. 


Mp(MoV/c 2 ) 


r^(McV/o 2 ) 




1 

2 
3 


778 ±7 
777 ±7 
773 ±7 


148 ± 14 
139 ± 13 
127 ± 13 


//// p = 0.47±0.07±0.12GeV 
|B/A|=0.81±0.08±0.20GeV~ 1/2 
ra = 5.7±0.4±1.5 



TABLE I: Parameters for different mass parameterizations. 



Cross Section 


STAR (mb) 


Ref [ 


] (mb) 


a 9 

^ xn,xn 
p 

a in.ln 


28.3±2.0±6.3 
2.8±0.5±0.7 


27 
2.6 


p(inc. overlap ) 
&x n,xn 
p 

^xn,On 

a 9 


39.7±2.8±9.7 

95±60±25 
370±170±80 




p 

total 


460±220±110 


350 



TABLE 11: Comparison to predictions from 
tainties are highly correlated. 



The uncer- 



parametrization [[22| . Here, T p = T ■ (M p /M^) ■ [(M ^ 
4m^)/(Mp— 4m^)] 3 / 2 is the momentum-dependent width, 
and f p is a fixed second order polynomial describing the 
residual background. The fit parameters are given in 
Table Q. The p° mass and width are consistent with ac- 
cepted values p3| ]; they were fixed to reduce the num- 
ber of degrees of freedom to obtain \B /A\, fi / f p , and 
ri. Our results are consistent with values found for the 
same parameterizations in jp — > p°p photo-production 
data(2l[|4). 

For coherent p° production accompanied by mutual 
nuclear break-up (xn,xn), we measure a cross section of 
<j(AuAu -> Au* xn Au* xn pV) = 28.3 ±2.0 ± 6.3 mb in the 
two-track event sample, by extrapolating the integral of 
the Breit-Wigner fit to full rapidity. By selecting single 
neutron signals in both ZDCs, we obtain a{ n i n /<J xn xn — 
0.097 ± 0.014, so a{AuAu -> Au* ln Au* ln p°) = 2.8±0.5± 
0.7 mb. Single neutron emission is predominantly due 
to Coulomb excitation and the subsequent decay of the 
giant dipole resonance. The ratio cr^„ i n /<jf^ nxn is con- 
sistent with crin,in/&xn,xn = 0.12 ± 0.01 found for mutual 
Coulomb dissociation at RHIC supporting that p° 
production and nuclear excitation are independent pro- 
cesses. 

At b ~ 2Ra coherent p° photo-production can overlap 
with grazing nuclear collisions, producing a low px p° 
accompanied by additional tracks. Additional tracks can 
also be produced at b > 2Ra from nuclear excitation 
by high energy photons. At present, we can not dif- 
ferentiate between these two processes. The coherent 
(xn,xn) p° sample increases by 40% when events with ad- 
ditional tracks are included. Accounting for this, we find 



a 



(inc. overlap) 



(AuAu -> Au* xn Au xn p°) = 39.7±2.8±9.7 mb. 
For (In, In) events, no additional p a candidates are found 
with higher track multiplicities. 

The major systematic uncertainties are in the 47r ex- 
trapolation (15%), acceptance and reconstruction effi- 



ciency (12%), luminosity determination (10%), and event 
selection (5%). The overlap region with grazing nu- 
clear collisions contributes 10%; it does not contribute to 
a in in) but a 10% uncertainty is due to the selection of 
the single neutrons. These contributions add in quadra- 
ture to 24% systematic uncertainty in the cross sections. 

The absolute efficiency of the year 2000 topology trig- 
ger is poorly known and does not allow a direct cross 
section measurement. From the two-track events, we ob- 
tain the cross section ratios a§, n xn /<To n 0n = 0.09 ±0.04 
and <jfj. n xn /cr xn On = 0.30±0.19. The uncertainties reflect 
the small number of (xn,xn) and (xn,0n) events in the 
topology trigger data. Grazing nuclear collisions do not 
contribute to a^ n Qn and a xn 0n , since they yield neutron 
signals in both ZDCs. From a(AuAu — > Au* xn Au* xn p°), 
we estimate <r{AuAu — > AuAup ) = 370 ±170 ± 80 mb, 
a(AuAu -> Au* xn Aup°) = 95 ±60 ± 25 mb, and the to- 
tal cross section for coherent p° production a (AuAu —> 
Au^AuWp ) = 460±220 ± 110 mb. Table || compares 
our results to the calculations of Ref. ||. The calcu- 
lation for a% n xn excludes grazing nuclear collisions; it 
is therefore compared to our value without the over- 
lap correction. Recent predictions [^5| are about 50% 
higher than in Ref. ^ without giving specific numbers 
for y/^w = 130 GeV. 

In summary, the first measurements of coherent p° pro- 
duction with and without accompanying nuclear excita- 
tion, AuAu — » Au*Au*p° and AuAu — > AuAup , con- 
firm the existence of vector meson production in ultra- 
peripheral heavy ion collisions. The p° are produced at 
small transverse momentum, showing the coherent cou- 
pling to both nuclei. The cross sections at ^/snn — 
130 GeV are in agreement with theoretical calculations 
based on the Weizsacker- Williams approach for large rel- 
ativistic charges, the extrapolation from p°-nucleon to p°- 
nucleus scattering, and the assumption that p° produc- 
tion and nuclear excitation arc independent processes. 
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